Heat capacity uncovers physics of a frustrated spin tube 



Nedko B. Ivanov^^^p Jiirgen Schnack^llj Roman Schnalle^, Johannes Richter^, Paul 
Kogerler"*, Graham N. Newton^, Leroy Cronin^, Yugo Oshima^, and Hiroyuki Nojiri^ 

^Institute of Solid State Physics, Bulgarian Academy of Sciences, Tzarigradsko chaussee 72, 1784 Sofia, Bulgaria 
^ Fakultdt fiir Physik, Universitdt Bielefeld, Postfach 100131, D-33501 Bielefeld, Germany 
^Institutfiir Theoretische Physik, Universitdt Magdeburg, P.O. Box ^120, D-39016 Magdeburg, Germany 
"^Institut fiir Anorganische Chemie, RWTH Aachen, Landoltweg 1, D-52074 Aachen, Germany 

^ Dept. of Chemistry, The University of Glasgow, Glasgow, G12 8QQ, UK and 
^Institute for Materials Research, Tohoku University, Katahira 2-1-1, Sendai 980-8577, Japan 

(Dated: July 11, 2011) 

We report on refined experimental results concerning the low-temperature specific heat of the 
frustrated spin tube material [(CuCl2tachH)3Cl]Cl2. This substance turns out to be an unusually 
perfect spin tube system which allows to study the physics of quasi-one dimensional antiferromag- 
netic structures in rather general terms. An analysis of the specific heat data demonstrates that at 
low enough temperatures the system exhibits a Tomonaga-Luttinger liquid behavior corresponding 
to an effective spin-3/2 antiferromagnetic Heisenberg chain with short-range exchange interactions. 
On the other hand, at somewhat elevated temperatures the composite spin structure of the chain is 
revealed through a Schottky-type peak in the specific heat located around 2 K. We argue that the 
dominating contribution to the peak originates from gapped magnon-type excitations related to the 
internal degrees of freedom of the rung spins. 
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Introduction-Spin tubes constitute a special class of 
quasi-one dimensional spin ladder systems characterized 
by periodic boundary conditions in the rung direction 
^U^. The magnetic compound [(CuCl2tachH)3Cl]Cl2 
is a geometrically frustrated triangular spin tube, the 
frustration being related both to the triangular arrange- 
ment of its rungs and to the twisted geometry of the legs, 
compare Fig.[lja). The relatively simple exchange path- 
way structure, described only by two dominant Heisen- 
berg exchange couplings IJ, [l5| , as well as the extremely 
weak exchange interactions between neighboring tubes 
renders [(CuCl2tachH)3Cl]Cl2 an excellent real mate- 
rial to study general properties of the spin-tube systems 
[1, 0, [13, 113, 13l ■ An appropriate spin Hamiltonian de- 
scribing the magnetic properties of this material reads 
as 



— ^ ^ Cn.Q • fn.Q-l-l (1) 
n— 1 a— 1 

where <t„c (a = 1,2,3) are spin-1/2 operators defined 
on the vertices of the elementary triangle denoted by in- 
dex n {n — 1, . . . , L). As depicted in Fig. [T] (b), the 
twisted spin tube may also be thought of as a three-leg 
ladder with periodic boundary conditions in the rung di- 
rection, where the parameters Ji and J2 are strengths 
of the rung and crossing (diagonal) exchange bonds, re- 
spectively. [(CuCl2tachH)3Cl]Cl2 is characterized by the 
parameters Ji/ks — 1.8 K and J2/kB — 3.9 K [isf . 
whereas the leg exchange constant J2 > 0, introduced 
for the sake of clarity in Fig. [1] (b), seems to vanish in 



this material. Figured] (b) clearly reveals the translation 
symmetry with one triangle per unit cell. 




FIG. 1: (a) Sketch of the twisted spin-tube system, (b) An 
equivalent spin model obtained through inversion of every sec- 
ond triangle. The nth elementary cell contains the spin-1/2 
operators fT„,i,o-„,2, and cr„,3. 

Two extreme scenarios for the ground state of Eq. ([T]) 
with antiferro mag netic couplings (Ji, J2 > 0) were out- 
lined in Ref. jlO|. In the case of dominating Ji cou- 



plings, the system effectively maps onto an effective spin- 
chirality model, where the additional chirality degrees of 
freedom appear as a result of the ground-state degener- 
acy of each individual triangle. On the other hand, for 
dominating J2 couplings the system maps onto an effec- 
tive spin-3/2 antiferromagnetic Heisenberg chain (AHC) 
containing additional biquadratic exchange couplings. In 
this Letter, we demonstrate by means of refined specific 
heat measurements, that the low-temperature proper- 



2 



ties of the spin tube material [(CuCl2tachH)3Cl]Cl2 re- 
produce the behavior of a spin-3/2 AHC characterized 
by the effective short-range exchange couphng constant 
Jeff = 2J2/3. Since the exchange interactions between 
different tubes are extremely small, the discussed com- 
^und also provides a rare example of spin-3/2 AHC [13- 
2l|. The experimental observables strongly suggest that 
if the system ever orders this should be much below 0.1 K. 
At elevated temperatures the measured specific heat ex- 
hibits a big Schottky-type peak located around T « 2 K. 
A detailed analysis - combining the semiclassical spin- 
wave approach with a number of numerical techniques 
such as the Quantum Monte Carlo (QMC) method, the 
Lanczos exact numerical diagonalization (ED), and the 
complete exact diagonalization [i^, [l^ - implies that 
the main contribution to the specific-heat peak stems 
from the lowest-lying gapped magnon excitations result- 
ing from the internal degrees of freedom of the composite 
rung spins. 

For the following discussions it is instructive to rewrite 
Eq. ([T]) in the form 



+ J2 Sr, 



V 



(2) 



Here s„ = (T„ 1 + (t„ 2 + o'n,3 is the rung spin operator 
related to the nth triangle. The number of Cu sites is 
— 3L, where L is the number of rungs. The interaction 
term V reads as 



V — (J2 ^ J2) 0-n,a ■ Crn+l,a , 



(3) 



tube material [(CuCl2tachH)3Cl]Cl2 are dominated by 
three types of low-lying spin excitations. Apart from the 
standard gapless excitations (characteristic of any half- 
integer AHC with short-range exchange couplings), im- 
portant contributions to the specific heat appear from 
two additional branches of low-lying gapped magnon- 
type modes, which are related to the chirality degrees 
of freedom of the local triangles. 

A qualitative picture of the low-lying spin excitations 
can be obtained already in the framework of the semiclas- 
sical spin- wave approach starting from the classical Neel 
configuration \St, —St, ■ ■ ■) , where St is the maximal value 
of the ^-component of the rung spin {St = 3/2 in the 
present case). Since the elementary cell contains three 
spin-S" variables, there appear three different branches of 
spin-wave modes 

E^{k,) - ASJ2 X (4) 
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— - [cos kx+cos{kx + ky)] 



where a — J1/J2, ky — 27rm/3 (m = 0,1,2), and S is 
the spin of a single site {S ~ 1/2 in the present case). As 
may be expected, the energy of the m = branch does 
not depend on the parameter Ji , since it is related to the 
dynamics of the cell spins s„ as a whole: 



Eo{k,j;) = Ws| sinfc^l, 



vs = (65)^ . 



(5) 



The above expressions reproduce the well-known semi- 
classical results for the dispersion relation and the re- 
lated spin- wave velocity Vg of a spin- (3 • S) AHC with the 
effective exchange constant J^tf — 2J2/3. 



where we have included an additional leg exchange in- 
teraction J2, see Fig. [1] (b). The last two equations ex- 
plicitly show that the condition J2 = J2 defines a spe- 
cial symmetric line in the parameter space on which the 
rung spin operators s„ are conserved quantities. On this 
line, the Hilbert space is decomposed into the invariant 
subspaces (sectors) [si, S2, . . . , sl], where the local spin 
quantum numbers s„ — 1/2, 3/2 are defined, as usual, by 
the relations = s„(s„ -1-1) (n = 1, . . . , L). In partic- 
ular, if the Hamiltonian ground state lies in the sector 
[3/2, 3/2, . . . , 3/2], then Eq. © for J2 = j'2 will describe 
a spin-3/2 AHC characterized by the effective exchange 
constant JeS = J2- Up to first order in the (formally) 
small parameter ( — J2)/ Ji , the effect of the interaction 
term V in Eq. ([2]) is reduced to a simple renormalization 
of the exchange constant in the effective spin-3/2 AHC: 
Jeff — i- Jeff = J2 + (J2 ~ J2)/3. At the special point 
J2 = 0, the result Jeff = 2J2/3 coincides with the first- 
order result of Ref. [l^l, which is obtained by another 
perturbation scheme starting from the limit | Ji | 3> J2 
(Ji < 0). 

Low-lying spin excitations- As argued below, the con- 
tributions to the low-temperature specific heat of the spin 
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FIG. 2: Spin-wave excitation modes (solid curves, Eqs. ([5]) 
and ([6}) compared with the lowest triplet excitations in two 
periodic clusters (symbols). 

The dispersion relations of the ky = 27r/3 and ky = 
47r/3 excitations can be expressed in the following form 



i?i.2(fcx) = J^' + 4^'^l sin^ (^k, T y) ' 



(6) 
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where the excitation gap A = 4S'J2^(1 — 3a/4)2 — 1/4 
corresponds to the lowest-lying modes at wave vectors 
kx = tt/S and 2tt/3. 

The solid curves in Fig. [2] show the dispersion rela- 
tions of the discussed spin- wave excitations. On the other 
hand, the symbols depict the positions of the lowest-lying 
triplet states, as obtained by the ED method for periodic 
clusters containing L = 8 and 12 unit cells 0. Apart 
from the finite-size effects related to the ED results, it is 
clearly seen that both methods qualitatively yield sim- 
ilar results. As a matter of fact, the discussed spin- 
wave branches may also be considered as one-dimensional 
analogs of the three spin- wave branches in the triangular 
lattice antiferromagnet. In this respect, the lowest-lying 
excitations at fc^; = 0, 27r/3 and 47r/3 in the spin tube are 
one-dimensional analogs of the three Goldstone modes in 
the triangular lattice antiferromagnet. 

More accurate estimates for the parameters of the ex- 
citation spectrum Vg and A can be obtained from an 
extrapolation of the ED results for L = 6, 8, 10, and 12 



unit cells. Using the approach of Ref. [18|, one finds the 
following estimates from the extrapolations of Eq{2tt/L) 
vs. sin(27r/L)/L and Ei{2tt/3) vs. 1/L (see Fig.^: 



vJkB = 10.06K = 3.87 



2.h 

3kB 



A/kB = 5.32K. (7) 



Interestingly, the extrapolation result 3vs/{2J2) — 3.87 
exactly reproduces the density-matrix renormalization 
group estimate for the spin- wave velocity of the spin-3/2 
AHC characterized by the exchange constant 2J2/3 [iSj. 
As already discussed, the same effective exchange con- 
stant (Jeff = 2J2/3) arises both in the first-order result 
for the effective spin model and in the semiclassical spin- 
wave approach. Below we demonstrate numerically that 
the specific heat of a spin-3/2 AHC with the exchange 
constant 2J2/3 excellently reproduces the experimental 
results in the low-temperature region T < 0.5 K. 
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FIG. 3: Extrapolation of the ED results giving the parameters 
of the excitation spectrum Vs (l.h.s.) and A (r.h.s.). 

Low-temperature specific heat-The specific heat. 
Figs. |4] and [SJ was measured at the Institute for 
Materials Research (IMR) of Tohoku University us- 
ing polycrystalline samples of the spin tube material 
[(CuCl2tachH)3Cl]Cl2. The solid curve in Fig. H depicts 
the specific heat of a spin-3/2 periodic AHC composed 



of 100 spin sites. The specific heat is evaluated by means 
of a QMC method employing the ALPS code [2J]. As 
an effective exchange parameter we used Jos = 2J2/3. 
As can be seen in Fig. 01 the QMC result reproduces the 
experimental data very well in the region T < 0.5. This 
means that even for = 0, when the parameter ( J2~^2) 
is definitely not small, the low-energy physics of Hamil- 
tonian H]) is described by the spin-3/2 AHC. Therefore, 
it turns out that the theoretically predicted biquadratic 
exchange term [J [l3| plays no role in the experimentally 
interesting region of the phase diagram characterized by 
the dimensionless parameter J1/J2 = 0.46. 



0.15 



^ 0.10 



0.05 



0.00 




O.i 0.6 

T (K) 



FIG. 4: Specific heat (per Cu spin) of [(CuCl2tachH)3Cl]Cl2 . 
The symbols denote the experimental values. The solid curve 
is the QMC result for a spin-3/2 chain of length L = 100. The 
dashed line provides the linear specific heat corresponding to 
the universal Tomonaga-Luttinger liquid form presented by 
Eq. ((SJ, by using the extrapolation result Vs = 10.06 K. 

Turning to the extremely low-temperature regime, 
we expect the universal specific-heat behavior of the 
Tomonaga-Luttinger liquid 



g(T) 

NkB 



ncT 
9^ 



(8) 



where the constant c stands for the so-called topologi- 
cal charge (c = 1 for a Tomonaga-Luttinger liquid state) 
and Vs is the velocity of the gapless spin excitations. As 
clearly seen in Fig.|4l already for T < 0.5 K the measured 
specific heat coincides with the QMC results and nicely 
extrapolates towards the universal behavior repesented 
by Eq. ([5]). The latter observations strongly imply that 
the spin tube compound [(CuCl2tachH)3Cl]Cl2 is char- 
acterized by a gapless Tomonaga-Luttinger liquid ground 
state. In a recent report. Nuclear Magnetic Resonance 
measurements also indicate a gapless spin state in the 
same material based on estimates for the extremely low- 
temperature part of the magnetic susceptibility [25 1. 

To explain the experimental data at intermediate tem- 
peratures around 2 K, we use the established structure 
of the low-lying excitation spectrum, compare Fig. [2l As 
a rough approximation, one may take i5i,2(^a;) ~ A and 
use the well-known expression for the specific heat of a 
two-level system (with the assumption that the excited 
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level has twice the weight of the ground state, i.e., r = 2 
in the following expression). Since the exact density of 
states is not known in the thermodynamic limit, we use 
an overall parameter A in order to fix the height of the 
Schottky peak: 

C ^ r(A/r)^cxp(A/r) 
NkB [exp(A/r) + r]2 ' ^ ' 

Notice that the position of the peak does not depend on 
the value of A. The expression for C{T), Eq. with 
A/fcs — 5.32 K is plotted in Fig. [S] by a thick curve 
together with the experimental data. One observes that 
it reproduces very well not only the position of the peak 
[Tm « 2 K) but also the behavior of the specific heat 
for T < Tm down to T « 0.7 K, where the contribution 
from the gapless branch Eolkx) becomes important. In 
addition, in Fig. [5] we also show the specific heat that 
results from complete diagonalizations of a few finite-size 
periodic clusters. As easily seen, the overall agreement is 
good, in spite of the pronounced finite-size effects. 
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FIG. 5: Specific heat (per Cu spin) of [(CuCl2tachH)3Cl]Cl2 . 
The symbols denote the experimental values, whereas the 
solid curve depicts the two-level approximation. The broken 
curves denote the specific heat for three complete diagonal- 
izations for finite sizes. 

Turning to higher temperatures, a few words are neces- 
sary to expain the discrepancies between the experimen- 
tal and numerical results for C{T). Notice that already 
at T « 2 K the phonon contribution begins to domi- 
nate the specific heat. Since the detailed phonon spectral 
density is unknown, the raw experimental data for C(T) 
was - as usually - corrected by subtracting a reasonable 
Debye-like specific heat contribution |26j]. 

Condusion-We demonstrated that the low- 
temperature specific heat behavior of the spin tube 
compound [(CuCl2tachII)3Cl]Cl2 suggests a Tomonaga- 
Luttinger liquid type ground state for this material, 
corresponding to an effective spin-3/2 antiferromagnetic 
Heisenberg chain characterized by the short-ranged 
exchange-coupling constant Jcg — 2J2/3. On the other 
hand, we argued that the main contribution to the 
observed Schottky-type peak around T « 2 K comes 



from the lowest-lying gapped magnon-type excitations 
resulting from the internal degrees of freedom of the 
composite rung spins. 
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